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We study the effects of solute interactions on osmotic transport through pores. By extending 
single-file, single-species kinetic models to include entrance of solute into membrane pores, we model 
the statistical mechanics of competitive transport of two species across membrane pores. The results 
have direct applications to water transport across biomembrane pores and particle movement in 
zeolites, and can be extended to study ion channel transport. Reflection coefficients, the reduction 
of osmotic fluxes measured using different solutes, are computed in terms of the microscopic kinetic 
parameters. We find that a reduction in solvent flow due to solute-pore interactions can be modelled 
by a Langmuir adsorption isotherm. Osmosis experiments are discussed and proposed. Special cases 
and Onsager relations are presented in the Appendices. 

I. INTRODUCTION 

Recent X-rayEI and electron crystallographic studiesBi have yielded structures of integral membrane proteins such 
as K + ion and water channels to near atomic resolution. Since many biological transport chantusk have specificity 
in allowing specific molecules to permeate and mediate simultaneous flows of numerous speciesiM These new data 
may offer insights that may help correlate structure to function. The pores spanning biomembranes as well as those 
found in zeolites (important in numerous industrial processes such as hydrocarbon separations and catalytic agents) 
are typically few Angstroms in diameter, and contain very few molecules at any time. Therefore, we will study 
osmosis-driven particle transport across such pores using simple one-dimensional lattice models. 

The next section reviews the linear phenomenological expressionsa and their parameters, allowing for a second 
interacting, competing species. The empirical parameters arise from macroscopic considerations only and are not 
derived fxom any microscopic pore structure. Section III formally defines a two-species model similar to that used by 
Su et alxi and Wang et alH Here, we assume a microscopic structure consisting of single-file pares. Unlike the single 
species case, where the statistics of one-dimensional particle dynamics are well understood,B~hil if solutes are allowed 
to enter a pore, no general analytical solution exists for the current across channels of arbitrary length. However, 
similar assumptions implicit in modelling single-species, single-file transport are used here: All nonlinearities except 
for those associated with the internal particle dynamics are neglected. Although our analysis is semiquantitative, and 
neglects some precise details about the specific particle-particle and particle-pore interactions, it provides a physically 
and mathematically consistent microscopic picture of nonequilibrium flow through single-file, solute permeable pores. 

In the Solutions and Discussion section, we plot the behavior of steady state flows under a variety of experimentally 
motivated conditions. Certain macroscopic conjectures of the role of solute size on osmotic permeabilityHj are re- 
examined. For short pores that allow partial and total passage of solute, the effective reflection coefficients (defined in 
Section II) and counterflow of solute are also computed and plotted in a number of illustrative graphs. An extension 
to longer pores is derived using an equilibrium multi-site Langmuir adsorption modelJi3 

In the Conclusions, we summarize our analyses and discuss interpretations of osmosis measurements. We argue 
that models which intuitively incorporate solute-membrane interactions into linear coefficientscJ are inconsistent with 
a virial expansion in the solute-pore interactions. Possible experiments and correlations among the rate parameters 
are proposed. For completeness, the lattice statistics, rate equations for two-sectioned pores, Onsagers relations, and 
analytic solutions in a special case are treated in detail in the Appendices. 
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II. PHENOMENOLOGICAL EQUATIONS 




FIG. 1. Schematic of osmosis and pressure driven flow through a membrane pore system. The reservoirs (L) and (R) are 
assumed infinite. The coefficients in the inset {£} = (a, /3, 7, 5,p, q), and {£'} = («' , f3' ,7' ,5' ,p' ,q') are microscopic conditional 
transport rates for solvent and solute respectively and are defined by the figure. The solvent(solute) mole fractions in the (L) 
and (R) reservoirs are Xl(x'l) an d Xr(x'r) respectively, and Xl,r + Xl r ~ 1- in a typical osmosis experiment, the membrane 
is impermeable to solute, {e.g. 7' or 8' — 0), and the hydrostatic pressure difference AP = Pr — Pl — 0. 



Wc first review the linear phenomenological equations describing two species transport across a membrane separating 
two infinite reservoirs (Fig. [l]). When a membrane is permeable to solute as well as solvent, their fluxes, J' and J, 
are coupled. Following Katchalsky and CurrarJj and defining the fluxes in terms of a volume flow Jy = vJ + v' J' 
conjugate to hydraulic pressure driving forces, and a relative flow Jo = J' /n' L — J/ul conjugate to diffusive driving 
forces, a set of linear Onsager relations are derived: 

J v = LpAP + L PD An 

(2.1) 

J D =L DP AP + L D AIl. 

Here, v(v') is the molecular volume of solvent (solute) in solution, rii(n^) is the solvent (solute) particle number 
density in reservoir (L). The osmotic pressure AIT = nxksTAx 1 , where tit is the total particle number density 
which we approximate as equal in the two reservoirs. The condition for zero volume flow (Jy = 0) requires AP = 
— (Lpd / L#)AI1. Thus, equilibrium across solute-impermeable membranes requires AP = All, and the reflection 
coefficients E = —Lpo/Lp = 1. For permeable solutes, the condition for zero volume flow yields a reflection 
coefficient S < 1 implying that a smaller hydraulic pressure is required to balance volume flow under a steady state 
osmotic pressure All, 

Jv = L P (AP - SAn). (2.2) 

Since Lpp> = —Lp (S = 1) for an impermeable solute, E can be measured by comparing Jy induced by permeable 
solutes to the maximum (Jy = vJ max ) induced by impermeable solutes. However, the experimental measurement of 

E = -^ = ^ + f},J/ (2.3) 
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(under isobaric AP = conditions) does not rely on the linearity assumed by (2.1), and nonlinear mechanisms may 
manifest themselves in E. Furthermore, as we will see, the current derived from an impermeable solute can also 
depend (to higher order in solute concentration difference A%') on how the solute particles partly enter and bind to 
the pore interiors. Thus, measurements of E may have more intricate dependences on solute-pore structure than is 
sometimes taken into account. 



III. TWO SPECIES ID EXCLUSION MODEL 

We now study flow mediated by one-dimensional channels. Single species exclusion modelsBlii~0 are extended to 
allow partial or complete entrance of solute particles inside a pore. This approach qualitatively models noncylindrical 
shapesEZl often found in biological pores and ion channels, i.e. flares, conical sectionsEJa, and vestibulescl in the pore 
structure. These wider sections may allow solute entrance and binding, and may render the pore permeable to solute. 
Ion specific channels are also believed to have a narrower section near the midplane which acts like a selectivity filter 
for only those ions which can bind or pass through. 

The one-dimensional chainjshawn in Fig. 1 which allows two types of particles {e.g. solvent and solute, or A and B) 
to enter and occupy any site.t3'E-3 The kinetic rates pictured in Fig. 1 are defined as follows. A solvent particle enters 
site i = 1 from (L) with probability per unit time dt or rate a\L if and only if site i = 1 is empty. This entrance rate 
is given by an intrinsic rate a times the mole fraction of solute particles in (L). Similarly, the entrance rate from (P) 
into the empty site i — N is given by Sxr- The exit rates from the pore to the reservoirs at site i = X,N provided 
they are occupied are 7, (3 respectively. In the pore interior, a solvent particle hops to the left with probability q only 
if the site i — I to the left is empty, transitions to an empty site to the right are denoted by p. If the particles do not 
experience pondermotive forces (electrostatic or gravitational forces) and are being transported across microscopic, 
uniform pores, p ~ q. However,, in general, the pores need not be uniform, and pondermotive forces (such as electric 
fields acting on charged particles) may also exist such that p / q. Completely analogous rates are defined for the 
second type of particle, or solutes, by primed quantities. The chain length is L = N£, where t is set by the size of 
the larger of the two types of particles. If a and a' are the approximate molecular sizes of the solvent and solute 
molecules, a > a'/2, then each site can contain at most one particle of either type. 



By considering transitions among all possible particle configurations, the kinetic steps defined by (Al) and (A2) in 
Appendix A can be written in terms of dynamical rate equations. When solutes cannot enter the pore (a' = 5' = 0), 
the solvent current J m ax a cross symmetric pores (p = q) in the presence of noninteracting solutes (i.e. a' = j3' = 7' = 
d'=p' = q' = 0) is exactly^ 

J ( N ) = p(«Xl/3 - l&XR) (3 n 

max[ ' (N-l)(axL + 7)(P + SxR)+p(axL + P + l + Sxny K } 

When the solute particles can enter the pore, no analytic expressions exist for J, J'. However, for illustration, 
consider a membrane so thin that N = 1, and only boundary entrance and exit rates are relevant. The single site 
within the membrane can be in only one of three possible states: empty (Pi), solvent filled (P2), and solute filled 
(P3). The probability fluxes are determined by the rate equations 

d t Pi = - (a XL + a'x' L + S X r + 5'x'r) Pi + (P + 7)^2 + (/?' + i)P* 

d t P2 =(axL+8 X R)Pi-(P + l)P2 (3.2) 

d t P 3 = (a'x' L + S' X 'r)Pi ~ (/?' + 7 0^3 

Upon imposing steady state (dtPi = 0), and normalization (Pi +P2 +P3 = 1), the averaged occupations (cr), (a 1 ) can 
be found and used in J(N = 1) = aXi(l — (f) — (cr')) — 7(0-) to find the steady state solvent current 

J(N = l) = (.aXLP-lSxR.)((3'+Y) f , 

1 ' ^Xl+P + 1 + 8 X r){P' +i) + {a'x' L + 5'xW + iy { '> 

The solute current J' is given by interchanging {£} = (a, j3, 7, 6) <-> {£'} = (a', /3', 7', 5') and xl.r Xl r m (HI)- 
Note that J(N = 1) has the same form as J ma x{N — I) appropriate for single-species transport except for an 
additional term in the denominator, 

(a'x'L +S'x'r)(I3 + i) 

/3' + 7' ' { ' 
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representing interference from states with solute occupation, a' = 1, (P3). The reflection coefficient defined for a 
single-site pore is thus 



S(JV = 1) 



1 + 



(aXL+l3+l+Sx R )(P' +1') 



(a'x' L +S'x' R )(P+l) 



1 



+ 



v J'(N = 1) 

i 1 Jmax{N = 1) 



(3.5) 



and represents two factors; reduction of solvent flow (first term), and solute backflow (second term). 

The microscopic activation energies E^(E'^), and the rate parameters {£} = (a,0,~/,8,p,q) and {£'} = 
(a',/3',7', 5',p', q'), are governed by the microscopic membrane-solvent (solute) interactions, and are approximated 
by the Arrenhius forms £ ~ £0 exp(— E^/ksT) and £' ~ £q exp(— E^/hsT). 

Kinetic equations can be readily solved for slightly longer pores. The states and transition probabilities for a two 
species, two section (N = 2) pore are enumerated in Appendix B. Currents across longer pores become increasingly 
tedious to compute as the matrix size increases as 3^ x 3^. The general transition rates for an N = 3 (27 x 27 matrix) 
pore are shown in Fig. @ and will serve as the basis for most our subsequent analyses. 



FIG. 2. Schematic of a three sectioned pore displaying all microscopic energies and transition rates. Separate energy land- 
scapes and occur for solvent(open circles) and solutes (filled circles). Exchange rates associated with wider pores are designated 
si,2. We take si = S2 = s for the same reasons q = p. 



Various effects of physical and chemical parameters, such as pore length and solute interaction dependences, are 
considered. We do not attempt to assign precise values to {£}, {£'}, but take the simplest, physically reasonable values 
to illustrate the relevant physical mechanisms. 



First consider solutes only in the (i?)-reservoir that can interact (i.e. partially enter) with a single-site (N = I) 
pore, but that encounters an infinite barrier while on the left side of the site (e.g. 5' 7^ 0,7' = x'l ~ 0)- Assume for 
simplicity a pore that remains microscopically symmetric with respect interactions with the solvent. A reduction of 
solvent flux manifests itself in the definition of S even though the pore is impermeable to solute and J' = 0. Thus 
E(N = 1, J' = 0) reduces to 






IV. SOLUTIONS AND DISCUSSION 



A. Solute Interacting, Impermeable Pores 



Xeff(N = 1) 



J(N = 1) 



Jmax(N = 1) 



(4.1) 




-1 -1 



= 1 + 



(2(a + 1) - a X ' R )(3' 



< 1 
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where a= a//3 defines an effective solvent-pore binding affinity.tS We find the reduction in osmotic solvent flux for 
solutes that enter the site, represented by £ eJ gs depends mainly on the solute equilibrium constant 5'x'r/P', and 
(a + Since occupation of the site by solute is required to influence (block) solvent flow, the S'x'r/P' behavior is 

not surprising. A low solvent affinity a further allows the site to be more likely occupied by solute. These competitive 
factors determine the likelihood solute enters the pore and hinders solvent flux. Note that S eJ y(x' fl — ► 0) — ► 1. 

The dependence on solute concentration x'r is not unexpected. Qualitatively, one power of x'r contributes the first 
order driving force for osmosis (the entropy of mixing); solute-solute (which we do not consider) and solute- membrane 
interactions must pome at higher orders of A\' via a virial expansion. Theories that suggest partial solute entrance 
into the porestSlij as a mechanism for reducing solvent flows require the presence of solute in the vicinity (adsorbed 
near) the pore mouth. The likelihood of this configuration will depend on the bulk solute concentration. Thus, these 
types of nonlinearities can be effectively built into a concentration dependent P os — P os (Ax').E3 In the Ax' — ► limit, 
P os approaches a constant (that depends only on how solvent travels through the membrane pore). The lowest order 
term is inde pend ent of solute identity. I— ■ 

Equation ( |4.l| ) arises from a Langmuir adsorption isothermlia determining the statistical fraction of solute- free (thus 
solvent conducting) states: 



J eff 



1-P eq {a' = l] 



1 + cxxl/P 



(4.2) 



where Z oc 1 + &Xl / P + $'x'r/0' is the partition function inc orporating the thr ee d istinguishable pore occupancy 
states of the entire pore-(R) reservoir ensemble. Equation (4.1) is recovered from (4.2). 

While a single occupancy pore may not be molecularly realistic, many of its qualitative physical characteristics 
remain relevant for longer, multi-occupation pores. We explicitly compute flows in two (Appendix B) and three 
sectioned pores by calculating the coupled 9x9 and 27 x 27 rate equations. For these longer pores, the reduction 
of Jy < Jmax due to decreased solvent flux will also depend on where along the pore the membrane first becomes 
impenetrable to solute. 
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FIG. 3. (a). Effective reflection coefficients E ej gr = J(N = 3)/ J ma x(N — 3), as a function of solute-pore attraction /3' , for 
solute impermeable pores (7' = 0, solid curves). The choice a = S = 8' — 0.01, j3 = 0.003 corresponds to a slightly solvent 
attracting pore. The three sectioned pore allows solute entry into N — 1,2, and 3 of its sections, schematically shown in (b). 
Also shown by the broken curves are £ Ej ff (this does not include the volume transfer by backflow of solute, treated in Section 
III) when the s olute is allowed to slowly permeate the pore, (f'//3' = 0.1,0.2). The curve labelled N = 10 is the equilibrium 
approximation (4.3) for a pore that allows solute entry into ~ 10 sites. Note the accuracy of (1.3) for N = 3, which is nearly 
indistinguishable from the exact result, (b) Schematic of asymmetric pores allowing varying degrees of solute entrance. 



Although solvent flux further decreases as solute is allowed to permeate deeper into the pore (allowing more solvent 
blocking configurations) solute-pore binding energies influence solute occupation exponentially. The sensitivity to 
solute-pore binding energy for various iV is demonstrated by Fig. ||](a), which shows the effective reflection coefficients 
for a three sectioned pore, T, e jj(N — 3, N a ) — J(N, N a )/ J max (N), as a function of the relative pore-solvent and pore- 
solute affinities. For concreteness we take dt = lOfs, such that (vt/£) ~ lps -1 sets p$ = q\ ~ 0.01 <C 1. Across pore 
sections that allow solute passage, we also assume p' = q' = 0.01. These are reasonable estimates when considering 
molecular length, mass, and energy scales under ambient conditions. Sections that do not allow solute passage are 
defined by p' = or q' = 0. The parameters used are s = a'x' L = 7' = 0, a = S = 8' = 0.01, and (3 = 0.003. These 
values correspond to a slightly solvent-attracting pore (a ~ 3.33), with identical intrinsic entrance rates for solvent or 
solute particles. The concentration x'r = 0.1/55.556 corresponds to a lOOmM aqueous solution in (R). The measure 
of relative affinity log 10 (/3//3') is varied by tuning the solute-pore attraction /3' from a' Q = 0.01 — > 0. N a < N is the 
number of sections into the pore the solute can enter from (R), as shown in Fig. |3|(b). For example, when N Q = 2, 
Pi = P2 = Qi = 0-01 an d q'i — 0. The upper three solid curves in Fig. [| represent the flux reduction from J m ax(N = 3) 
due to solute entrance into one, two, and all three sites, corresponding to q[ = 0, q% = 0, and 7' = respectively. 

The decrease of E eJ y for larger iV (the number of sites open or accessible to solute) can be physically understood 
in terms of a multisite Langmuir adsorption isotherm similar to that applied to the one-site model. The probability 
that there are no solute particles in any of the N a accessible sites is 



Zeff(N,N ) 



1 + Sxr/P 



i + Sxr/P + S'x'r/P'. 



N„ 



(4.3) 



This Increasing N , enhances the probability of nonconducting pores and can also be thought of as an additional 
entropic factor fesT In N a in free energy difference favoring nonconducting (solute adsorbed pores) over conducting 
states. Equation ( ll.3 j ) is valid provided J is too small to significantly affect equilibrium occupancies. For the parameters 
used, the e stim ate (pl.3|) for S(A r Q = 3) is nearly indistinguishable from the exact solution. The effective reflection 
coefficient (4.3) for N > N a = 10 (7' = 0, s = 0) is also indicated in Fig. 0. Note that the estimated relative flow 



reduction, ( |4.3D , is independent of the total pore length N. 
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FIG. 4. Solute concentration dependences of Y, e g are strongest highly attracting solute-pore species (small 0') . The N 
dependence is negligible with under these parameters, consistent with the multisite Langmuir adsorption model (4.3). 



Figure |] shows the dependence of £ e jp (iV, N a = 1) on solute concentration for N = 1, 2, 3 at various solute bindings 
(3' = 0.00003, 0.0003, 0.003. All other parameters are identical to those used in Fig. ||. The reflection coefficient 
decreases for larger solute concentration and smaller j3' as expected. The independence of (4.3) on N is also confirmed 
by the curves in Fig. ^, which are each three indistinguishable curves corresponding to N = 1, 2, and 3. 



B. Solute Permeable Pores 

When solute can steadily pass from (R) to (L) on time scales of the measurement ,0 the solute flow can nonnegligibly 
contribute to total volume flow. This flow comprises the solvent flu x te nding in one direction and solute flux tending 
in another. Therefore, the measured reflection coefficient S (Eqn. (2.3)) can even be negative. The solvent current J 
can also be independently negative, i.e. driven back by a strong counterflowing solute. Negative solvent flow is shown 
by the broken curv es in Fig. where the pore is permeable to solute (7' > 0). The reflecti on c oefficients measured and 
defined by S eJ y or (2J3) can thus be negative. However, for the TV = 1 case, we see from (3J3) that increasing 7' from 
zero can actually increase solvent flow J. This occurs particularly as x'l an d is the consequence of the additional 
route for emptying the solute from the membrane site to (L), increasing the probability for a solute- free conducting 
state. This "backside" exiting effect will only occur for very short pores where a single kinetic step governed by 7' 
(in addition to /?') renders the pore conducting and will disappear as N increases. 

The total volume flux, Jy = v J + v'J' for N = 1 found from ( |3.3| ) and J' can also increase as 7' increases provided 
(dJv/d-f 1 ) > 0, which yields (for x'l = = <^ an d (3 = 7) 



I 
0' 



> 



2(a + 13)- a X ' R 
a X ' R (v/v>) - 28'x' R 



(4.4) 



Thus, a single site pore that is impermeable (7' = 0) but allows solute entry only (5' ^ 0) can actually increase net 
volume flow as it is made slightly permeable (7' > 0). This anomalous behavior is more prevalent when v' /v and X ' L 
are small, since the volume backflow due to solute would be small and solute re-entrance rates from (L) vanish. 
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FIG. 5. Currents for various internal solute barrier heights occuring at different positions. iV = 1,2,3 correspond to 1,2,3 
easily solute-accessible sites within an iV = 3 pore. The quantity x\ = q[, q' 2 , 7' measures the hopping rate to the left past this 
site. 
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In Fig. [| we relax the impermeable solute constraint and allow for passage of solute from (R) 
the rate limiting internal barrier is varied relative to (3' at various internal positions i, where x\ — q\ , q' 2 , or 7 
X 'i/P' ~~ * limit of J and J' approach those expected from simple osmosis resulting from an impermeable (although 
pore-entering) solute, as studied in Fig. 0. As impermeability is relaxed, J decreases, while J' becomes negative. 
The effects are most pronounced for pores that have the largest number of easy entrance sites (e.g. N Q = 3). In fact, 
for solute-binding pores where q[/ [3' > 1, and larger N, the negative solute flux can actually drag, or pump solvent 
backwards (N — 3 curves in (b) 



for 7 7/3' > 1). 
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FIG. 6. (a). Solvent and solute fluxes as functions of relative pore-solute, pore-solvent binding affinity log(a'/3/a/3') (with 
a,/3 fixed). To the left of the discontinuity, (defined by a' — 0.003, (3' = 0.01 and most evident for N — 1) the pore is 
solute-repelling, and a' is being varied, while to the right, the pore is being made increasingly solute-attractive by decreasing 
. (b). The associated total pore occupancies. 



Negative solvent flux is shown in Fig. Here, we assume a molecularly symmetric (with respect to both solvent 
and solute) pore such that j3' = 7', (3 = 7 and q = p is constant along the pore, as assumed in the specific dynamical 
model presented in Appendix A. Figure 0(a) shows the currents J(N — 3),J'(N = 3) as functions of the solute-pore 
affinity a' / j3' expressed in terms of the logarithm of the relative (to solvent-pore binding fixed at a = 0.01, (3 = 0.003) 
affinities. The discontinuity at \og(a' /(3' Q ) — log(0. 003/0. 01) indicates the crossover from solute-repelling to solute- 
attracting pores. Only at intermediate solute affinity, and longer pores (N > 3), does solute current J' < sweep 
solvent in a direction opposite from that expected from simple osmosis. Pores that repel solute rarely contain such 
particles that drive solvent up its chemical potential gradient A/x. However, pores that bind strongly to solute are 
choked off, and both J and J' diminish. Fig. ^(b) shows the associated total pore occupancies (<Ji) and ( a i)- 
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FIG. 7. Solvent and solute fluxes across symmetric solute permeable pores as function of solute concentrations, (a) \'l = 
as Ax' = x'r i s varied. The lower three thin curves and upper thick curves depict solute (J' < 0) and solvent flow respectively. 
The parameters used correspond to a solvent attracting pore with strong solute binding (solid curves), intermediate binding 
(dotted courves), and solute repelling (dashed curves) pores, (b) The currents as absolute concentration \'l m varied, with 
A X ' = 0.1/55.556 fixed. 



Figure 0(a) shows the dependences of J(N = 3) and J'(N = 3) on solute concentration difference A^' with \'l = 0- 
For strong solute-binding pores (low f3' , solid curves), the solvent current J(N — 3) < due to a strong solute 
backflow J'(N = 3) < 0. In Fig 0(b), we vary x't while keeping A%' = 0.1/55.556 fixed. The higher absolute solute 
concentrations in this case decrease solvent occupation allowing a larger current in the same direction as osmosis, and 
solute is swept in the same direction J' > 0. However, for too strong a solute binding, increasing solute occupation 
eventually clogs the pore, precluding efficient solvent flow (solid curves). Here, at low absolute concentrations, J,J'<0 
indicating that solute is pushing solvent back against the direction expected from simple osmosis. 
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FIG. 8. Currents across a symmetric, N — 3 pore as a function of si = S2 = s. Concentrations used are 
Xi =0,xb = 0.1/55.556. 

Effects of "slip" si = S2 = s (see Fig. |[) are explored in Fig. ^. The magnitudes of both J(N — 3) and 
J'(N — 3) are increased as s increases since the countcrflowing currents are occasionally allowed to more easily pass 
through each other. This variation is strongest for smaller /3' (strong solute-binding pore) since large solute occupation 
configurations would be most influenced by slip. Note the near symmetry between J and — J' when (3 — /3' = 0.003. 



V. SUMMARY AND CONCLUSIONS 



We have studied a commonly measured figure of merit, the reflection coefficient £ (which describes the reduction in 
fluid volume transport in the presence of permeable or membrane interacting solutes). The single site model of Kalko et 
al.t3 is extended to include longer pores and the effects of solute penetration and binding. The currents and reflection 
coefficints (Figs. |||^) were computed in terms of microscopic molecular parameters in some limiting, yet illustrative 
cases. Solutes that enter and bind to pore interiors block and reduce solvent flux without being completely permeable, 
yieldins-jw, effective reflection coefficient Tj e g < 1. Although higher order corrections to van't Hoff's law have been 
studiedtaS and is well understood in terms of virial expansions of solute-solute interactions in the reservoirs, solute 
fouling arises from higher order concentration terms resulting from solute- membrane/pore interactions. Although the 
phcnomenological Onsager coefficients Lp 1 Lpp> = Lp,p, hp, can be qualitatively approximated from our plots of J and 
J', care mus t be exercised when interpreting measurements since solute interference effects can be highly nonlinear, 
when Eqns. 



2.1 



do not apply. The magnitudes of the nonlinear effects depend on solute-pore interaction energies as 
well as the number of configurations that are solvent blocking, both of which can be inferred from microscopic models 
of the particular pore-solute system being considered. 

For a one-site pore that allows solute entrance but remains solute-impermeable "solute interference" nonetheless 

reduces the solvent flux J by a factor cx XR ex P (~E's< + E'^/kpT , where x'r 1S the solute mole fraction, and E' s 
and E'p are the solute-pore "in" and "out" energetic barriers. Because the nonlinearities are most pronounced for 
large S'/f3', e.g. a very large Eg or a highly solute-binding pore, strong solute binding leads to the most effective 
solvent flow reduction. Measurements of J vs. x'r ffl^h different solutes should have the same slope as x'r ~> 0. 
If different solutes can indeed yield different slopes,Ej the interactions are too strong for the linear limit to have 
been reached. Extreme cases are "solutes" such as mercurial compounds p-chloromecuribenzoate (PCMB) that can 
reversibly destroy aquaporin pore function at very low concentrationsEj due to strong binding; although here, the 
binding is so strong that pore structure is presumably altered. Exnerirnentally, the solute concentration dependence 
on reflection coefficients has been observed in various systems .EfO'LJ 

One possible experiment to probe solute blockage effects in asymmetric (with respect to solutes) pores and flow 
rectificationta is to compare the currents J and — J associated with solutes in the (L) and (R) reservoirs. For a 
symmetric pore, the magnitudes of the currents will be equal. However, for asymmetric pores, flow induced by solutes 
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in (L) will be different from that generated by solutes in (R), the difference being manifested in higher order terms 
in solute concentration. 

If the pore is permeable to solute (J' ^ 0), and E is defined by measuring the total volume flux Jy = v J + v'J', a 
negative reflection coefficient E < can even occur when the counterflowing solute adds to a possibly counterflowing 
solvent, or overcompensates for the co-flowing solvent. For solute-permeable pores, correlations between volume 
flow and solute size typically used in macroscopic descriptions, when extended to the microscopic regime considered 
here, are not straightforward since Jv depends directly on v' as well as on the rate parameters {£'}. The notion of 
molecular size alone, based on van der Waal's radii for example, neglects pore-particle attractions and nonadditive 
interactions among molecules such as hydrogen bonding. For example, urea may easily enter pores based on its small 
size (and thus have a larger N a ), but participates in strong H-bonding in aqueous solution, implying a very small 
a' / (3 1 (solute-repelling pore) since H-bonds need to be broken prior to pore entry. Urea might then be expected to 
yield a larger measured reflection coefficient than that expected from its molecular size aloneHD H-bonding in the bulk 
solution may also significantly change D' and affect the interpretation if unstirred layers are important. Furthermore, 
the additional solute interactions giving -cise to E eJ y,E provide additional energy scales which can yield apparently 
non-Arrhenius temperature dependences r3 

The treatment given suggests that correlations between thermodynamic solute-solvent parameters and pressure 
or osmosis driven flows can be partly experimentally determined. Specific volume changes of bulk solute-solvent 
mixing, and enthalpies of mixing determine v'/v and constrains {£,£'}■ For example, a solute with a negative 
enthalpy of mixing (with solvent), AH m , would have a smaller S' /[3' than a solute-solvent pair with larger or positive 
AH m . Equilibrium measurements of relative solvent and solute absorption in pores also give relative measures of 
the parameters <5'//3'.-and a/ (3. Further refinement of the parameters {£}, {£'} may be obtained with—more precise 
molecular mechanicsE3 or molecular dynamics dynamics (MD) using appropriate molecular force fields.E3 Direct MD 
simulation of osmosis across very thin membrane pores (single atomic layers) have even been performed for high solute 
concentrations^ Nonequilibrium MD simulations have also been used to compute dynamic separation factors for two 
species flows£i Monte-Carlo simulations may also reveal qualitatively interesting behavior in transport across longer 
channels. 
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APPENDIX A: ID CHAIN DYNAMICS 



For arbitrarily long pores, the dynamical rules of particle transport incorporating exclusion interactions can be 
formally implemented on a discrete lattice. If 0^(0^) <E {0, 1} denote the solvent (solute) particle occupation at site 
z, the instantaneous rules for the microscopic transport of the two species in the pore interior are mathematically 
defined by 

Ji,i+1 = M(l - - Oi+l) - - <7i)(l - CT-) 

(Al) 

J( i+1 = p'aKl - a i+1 )(l - a' t+1 ) - q'a' i+1 (l - a^l - oQ 

where (p,q,p' ,q') € (0, 1). The time variable has been suppressed for notational simplicity. These rules incorporate 
excluded volume into the dynamics of both solvent and solute. Equations (Al) are supplemented with the boundary 
transition rates 



(A2) 





= a(l - 


< t 1 )(1- 


- O - 701 




= (3cr N - 


-5(1- 


o^X 1 -Ojy) 


T< 


= a'{l- 




-<7j -701 


Jl 


= $'<t' n 




-Ojv)(l-0-^) 
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Upon taking the time or ensemble average of the above equations, and assuming the rate parameters are independent 
of the dynamical variables <t,<j', we formally obtain steady state flows that depend on {£,£'}, the dt — ► limit of 
the averages of {£,£'}• However, ,-tha, time averaged currents are not easily solved due to higher order correlations 
among a, a' such as ((jj+icr^), etc.t2Mi-3 Nonetheless, transport across very short (small N) pores can be solved using 
rate equations defined by the rates {£}, {£'} found from the mean of the distribution of {£}, {£'} E {0, 1}. 



APPENDIX B: TWO SPECIES TWO SITE MODEL 



Consider a membrane thick enough to accommodate a pore that can simultaneously fit two particles, solvent or 
solute, along its axis. We assume that the particle-particle interactions are dominated by excluded volume effects 
such that longer ranged microscopic interactions that cooperatively affect the particle hopping rates can be neglected. 
At any instant, either of the two partitions in the pore must be in one of three occupation states, empty, solvent (open 
circle), or solute(filled circle). Therefore, there are 3x3 = 9 microscopic states denoted by 
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(Bl) 



where the transitions among the probabilities are given by the coupled rate equations 

d t P l 



A l3 P J= Q 



(B2) 



with the transition matrix defined by A = 
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(B3) 



where q;t(<5t) = (xxl + (x'x'l(Sxr + S'x'r)- The coefficients are defined by the elementary steps shown in Figure 
|l|; a, P, 7, S and their primed analogues carry the same physical meaning as in the one site model. The quantities 
q(p),q'(p'), and s(s') define the microscopic rates of transitions P 2 — > ^3(^3 — > P2), P5 — > Pe{P& — ► P5), and 
P§ — > P-;{P-i — » Pg); respectively. Probability distributions Pj are found by solving Eq. (B2) with (B3) and the 
constraint Pj = 1. Form these, the steady state solvent and solute currents are 



J = (xxl{Pi + P 3 + Pa) - 7(P 2 + Pi + P 8 ) 

J' = «'x' L (P2 + P3 + Pe) - 7'(Ps + Pt + P 9 ). 
An analogous set of equations hold for the kinetics of a three section, two-species pore. 



(B4) 



APPENDIX C: N SITE, TWO-SPECIES CASE 



We generalize the single section model to N sections. Under special cases, analytic expressions can be found for 
general N. Using notation describing the one section, the time averaged fluxes of a and a' particles are 
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(CI) 



J = p{(j l - cxi+i) + (p- s){a l+1 a' z - (Xia' i+1 ) 
J' = p'ia'i - a' l+1 ) + (p' - s)(a l a' l+1 - cr l+ i^), 
since (cicr^ . . .) = 0. The steady state volume flux Jy = vJ + v'J' in the interior and at the boundaries are 
Jv,i = vp{(Ti - + v'p'ia't - a' i+1 ) + [(p - s)v - (p' - s)v'] (cr i+1 cr< - <Ti<j' i+1 ) 

Jv,i = (va XL + v'oc'x'l) - QM - Q'M (C2) 
Jv,n = -(vS X r + v'S'x'r) + Qr{ctn) + Q' R {v' N )- 

where 



(C3) 



The quadratic terms in Jy,i cancel provided 

(p - s )v = (p' - s)v'. (C4) 
Under this special condition, it is useful to define af = pvoi + p'v ' a\ so that the current can be summed as 

The volume flux Jy is completely determined when the boundary currents depend only upon a + , which is the case 
provided 

Ql = Qr = VP_ (CF , 
Q'l Qr W 
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Q'r 


= vSxr - 


-v'S'x'r^ 


-xfp 



The above relationships, along with (C4) provide constraints among a', 7', i5' given a, /3, 7, J. In the simple case 



of no-pass pores, s = 0, a simple scaling exists; 

p'=p^ 0=^ 7' = 74- (C7) 

where </> = v/v' . These special conditions imply larger (in the aqueous bulk phase) solute are more sluggish in hopping 
between sections and into and out of the pore, and are not qualitatively unreasonable. In this illustrative case, the 
volume flux is 

, , m -yjvS - i'6')Ax' + v(a/3 - 7 3)x' £ + v'ja'P - jSJXx rrw 
(N - 1)QlQr+Pv{Ql + Qr) 

an implicit function of AP. In the x'l = 0> isobaric ( AP oc a(3 — j5 — 0) limit the sign of J' is determined by vS — v'S' . 

APPENDIX D: ONSAGER RELATIONS FOR ONE SITE MODEL 

We have only explicitly considered osmotic flow between isobaric reservoirs (L) and (R). However, when a hydraulic 
pressure is applied, the enthalpies between the two baths will generally be unequal, and two independent thermody- 
namic driving forces may exist. For microscopically symmetric pores (a = S,a' = S'), the differences AP and AE' can 
be readily related to hydrostatic pressure differences. Expanding the solutions ( |3.3[ ) and J'(N = 1) about Ax' = 
and Pl = Pr, we find 
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J = P os Ml - L P AP 



J' = -P^An - L' p AP 



(Dl) 



where the coefficients in terms of the microscopic kinetic parameters are 

P, 



arfy' 



P' 



2nTkBT(aj'xL + ol'^x'l + 77') 
a'77' 

2n T k B T(aj'xL + a'lx'h + 77') 

«77 xlv 

2k B T(aj'xL + ol'^x'l + 77') 

a'liXL^' 



2k B T(a-y'xL + a'jx' L + 77') ' 
where we have assumed an attractive pore and linearized the exit rates using their Arrhenius forms according to 

and the Maxwell relation (dEp/dPp) = v. 



Upon forming the conjugate flows and using the definition (2.1), 

wYJ'XlV 2 + a'77'xiu' 2 



Lp = 



2k B T(aixh + a'jx' L + 77') 



«77 1; - a 77 w 

L PD - L DP - — — ■ — (U4j 

2n T k B T(aj'xL + a'~/x L + 77 J 

ajj'xL - ol'h'x'l 



2n 2 T k B T{aixL + a'lx' L + 77') 



which manifestly satisfies Onsager reciprocity. Corresponding expressions for longer pores are very cumbersome but 
can be evaluated with e.g. Mathematica. 
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